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ABSTRACT 

Galaxies in Hickson Compact Group 91 (HCG 91) were observed with the WiFeS 
integral field spectrograph as part of our ongoing campaign targeting the ionized gas 
physics and kinematics inside star forming members of compact groups. Here, we 
report the discovery of H II regions with abundance and kinematic offsets in the other¬ 
wise unremarkable star forming spiral HCG 91c. The optical emission line analysis of 
this galaxy reveals that at least three Hll regions harbor an oxygen abundance ^0.15 
dex lower than expected from their immediate surroundings and from the abundance 
gradient present in the inner regions of HCG 91c. The same star forming regions are 
also associated with a small kinematic offset in the form of a lag of 5-10 km s“^ with 
respect to the local circular rotation of the gas. Hi observations of HCG 91 from 
the Very Large Array and broadband optical images from Pan-STARRS suggest that 
HCG 91c is caught early in its interaction with the other members of HCG 91. We 
discuss different scenarios to explain the origin of the peculiar star forming regions 
detected with WiFeS, and show that evidence point towards infalling and collapsing 
extra-planar gas clouds at the disk-halo interface, possibly as a consequence of long- 
range gravitational perturbations of HCG 91c from the other group members. As such, 
HCG 91c provides evidence that some of the perturbations possibly associated with 
the early phase of galaxy evolution in compact groups impact the star forming disk 
locally, and on sub-kpc scales. 

Key words: galaxies: evolution, galaxies: individual (HCG 91c), galaxies: interac¬ 
tions, galaxies: ISM, ISM: Hll regions, ISM: abundances 


1 INTRODUCTION 

The environment is known to be a key factor influencing 
the pathways of galaxy evolution. However, the interconnec- 
tivity between the large scale, environmental mechanisms 
(e.g. gravitational interactions, ram pressure stripping) and 
the internal galactic processes (starburst, star formation 
quenching, nuclear activity, and both inflows and outflows) 
remains poorly understood. Part of the issue lies in the fact 
that galaxies are spatially extended objects, but are often 


* Email: frederic.vogt@anu.edu.au 


represented with a series of ensemble properties in single¬ 
spectrum studies. The advent of integral field spectroscopy 
(IFS) has opened a new way to study galaxies and their 
evolution on a spatially resolved basis at optical and near- 
IR wavelengths. The wealth of IFS surveys that have been 
undertaken in recent years provides direct evidence of the 
scientific potential of spatially resolved studies of a statis¬ 
tically significant number of galaxies. Such surveys include 


SAURON (Bacon et al. 

200] 

L de Zeeuw et al. 2002 Em- 

sellem et al.||2004|), PINGS 

(Rosales-Ortega et al. 20101, 

ATLAS"^^ (Cappellari et al.|2011|), CALIFA (Sanchez et al. 

2012), SAMI (Groom et al. 

2012), VENGA (Blanc et al. 
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2013) 


2014 


S7 ( Dopita et al. 
Bundy et al.||2014] 


^2014a|b | 


and MaNGA (Drory et al. 


In this article, we continue our analysis of the ion¬ 
ized gas physics and kinematics associated with star form¬ 
ing galaxies inside compact groups started in |Vogt et al.| 
(2013). Compact groups are isolated structures by defini¬ 


tion, and are as such often described as perfect laboratories 
to study the consequences of strong, multiple, simultane¬ 
ous gravitational interactions on galaxies. Compact groups 
are intrinsically less crowded than clusters (although they 
can be denser, see Hickson et al.|19^ ), so that the detailed 
structure of their large scale environment can in principle be 
better understood. In groups, galaxies are to a first approx¬ 


imation mostly subject to gravitational interactions (Coziol 
fc Plauchu-Frayn|2007 ), but the role of other processes such 


as ram-pressure stripping remains uncertain. 

Indeed, the presence of a hot halo inside some compact 
groups has been confirmed observationally with a Rosat sur¬ 
vey by Ponman et al.|(|1996 ). The analysis of Chandra ob¬ 
servations by Desjardins et al. (2013) showed that in some 


compact groups (and unlike in clusters), the diffuse X-ray 
emission is associated with individual galaxy members. Yet, 
other X-ray bright and massive systems have been found to 
match the X-ray scaling relations of clusters, and may be 


representative of an evolved state of compact groups (Des- 
jardins et al.|[20T4 ). Galaxies in compact groups are found 


to be HI dehcient (when compared to isolated galaxies with 
similar characteristics). [Verdes- Montenegro et al.| ( 200 1| re¬ 
ported a correlation between HI deficiency and X-ray emis¬ 
sion, suggesting galaxyohot IGM (Intergalactic Medium) 
interactions as a possible origin for the observed Hi defi¬ 
ciency. Yet, Rasmussen et al. (2008) report in a Chandra 


and XMM-Newton study the non-detection of X-ray emis¬ 
sion in 4 (out of 8) of the most Hl-deficient compact groups 
of Verdes-Montenegro et al. (2001), suggesting that galaxy 


ohot IGM interactions (e.g. ram pressure stripping) may 
in fact not be the mechanism driving the observed HI defi¬ 
ciency. The detection of warm H 2 emission at levels incon¬ 
sistent with X-ray heating or AGN activity in 32 out of 74 


galaxies in 23 compact groups has lead Gluver et al. (2013) 


to propose that shocks induced by galaxies interacting with 
a cold IGM may be present in many compact groups. The 
Stefan’s Quintet compact group and its intergalactic shock 


represents an extreme example of this mechanism ([Appleton 
et al.||2006 Gluver et al.|2010 ). 


Within the larger scheme of galaxy evolution, it has 
been suggested that galaxies first start to evolve in compact 
group-like environments, before further processing in clus¬ 


ters. This is usually refereed to as “pre-processing” (Gortese 
et al.||2006 Vijayaraghavan Ricker||2013 ). Pre-processing 


mechanisms may not be restricted to compact groups and 


could be active over a wider range of environments ( Gy- 
bulski et al.||2014' ), but compact groups in particular favor 


a more rapid evolution of galaxies compared to the field. 
The existence of a gap in the mid-infrared color distribu¬ 
tion for galaxies in compact groups compared to the field 
has been interpreted as the ability of these dense environ¬ 
ments to rapidly transform star forming galaxies into pas¬ 
sive ones (Johnson et al. |2007[ [Gallagher et al.|200^ j Walker | 


et al.|201Ql' 2012). Based on the UV properties of galaxies in 

groups, Rasmussen et al. (2012) report that at fixed stellar 


~40 per cent less than in the field; a trend best detected 
for galaxies with masses < 10^ Mq. In the case of elliptical 
galaxies, jde la Rosa et ah (2007) found their stellar popu¬ 
lation to be older by 1.6 Gyr and more metal-poor by 0.11 
dex in [Z/H] compared to similar galaxies in the field, which 
they interpreted as the signature of truncated star formation 
in these systems. Yet, most galaxies in compact groups ap¬ 


pear to follow the B-band Tully-Fischer relation (Mendes de 
Oliveira et al.|2003 Torres-Flores et al.|2010 ), as well as the 


K-band Tully-Fischer relation ( Torres-Flores et al.||2013 ). 

Whether quenched or enhanced, changes in the star 
formation activity in a galaxy are merely a consequence 
of the formation, destruction, processing, spatial redistri¬ 
bution and local density variations of the molecular gas (i.e. 
the fuel) within the system. Therefore, understanding the 
large scale gas flows in compact groups is key to a bet¬ 
ter understanding of the evolutionary pathways of galax¬ 
ies in these environments. For example, jVerdes-Montenegroj 
et al. (2001) proposed an evolutionary sequence for com¬ 


pact groups, where the Hi gas distribution is first associ¬ 
ated with the individual galaxies, before gravitational in¬ 
teractions strip it apart, resulting in no Hi gas left in the 
galaxies, or possibly (but less frequently) in a large common 
envelope. 

Here, we describe the use of the WiFeS integral field 
spectrograph to study the physics and kinematics of the ion¬ 
ized gas in star forming galaxies inside compact groups. This 
approach is complementary to group-wide Hi observations 
(e.g Verdes-Montenegro et al.||2001 ) in that we focus on a 
different phase of the interstellar medium (ISM), and zoom 
in on the stellar and ionized gaseous content of the galaxies 
themselves. By targeting galaxies at z<0.03 with a spectral 
resolution R=7000 and a spectral pixel (spaxel) size of 1 
arcsec^ (with natural seeing), we gain access to the physics 
of the ionized gas on scales of ~1 kpc and simultaneously 
resolve velocity dispersions down to cr > 20 km s“^: an ideal 
combination to detect and study localized consequences of 
the compact group environment on galaxies. 

Our targets are drawn from both the Hickson and 
Southern Gompact Groups. The Hickson Gompact Groups 
(HGGs, [Hickson 1982a|b ) represent a class of compact 
groups that was first identified on Palomar Sky Survey red 
prints as groups of four or more galaxies with specific bright¬ 
ness, compactness and isolation selection criteria. Most were 


later on confirmed to be gravitationally bound (Hickson 


et al. 


1992 


Ponman et al. 


1996). Along with their south¬ 


ern sky cousins (Southern Gompact Group or SGG, Iovino[ 


mass, the specific star formation rate of group members is 


2002), the proximity of these low redshift systems (median 
0.03, Hickson et al.[[T9^ ) make them optimum targets 
for spatially and spectrally detailed observations. 

In Vogt et al. (2013), we focused on NGG 838 in 
HGG 16, confirmed the presence of a large scale galactic 
wind and characterized it as young, photoionized and asym¬ 
metric, very much unlike the near-by shock-excited galactic 
wind of NGG 839 ( Rich et al.|2010 ). In this article, we turn 
our attention to HGG 91c, a star forming spiral of which 
the basic characteristics are given in Table In an effort to 
tie the large scale structure of the group to the local envi¬ 
ronment within the optical disk of HGG 91c, we supplement 
our optical IFS observations with HI observations from the 
Very Large Array ( VLA) and with broadband optical im¬ 
ages from the Pan-STARRS survey. Effectively, we build a 
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Abundance and kinematic anomalies in HCG 91c 3 


Table 1. Basic characteristics of HCG 91c 


-27°45'00" 


46'00" 


l'=30kpc 


HCG 91c 


47'00" 


^ 48'00" 


49'00" 


50'00" 


51'00" 





09 ^ 00 ^ 22 '’ 08 "’ 54 ^ 


Property 

Value 

Ref. 

Names 

HCG 91c 



ESO 467 - G 013 


R.A. [J2000] 

22^09"'07.7^ 


Dec. [J2000] 

-27°48'34" 


Redshift 

0.024414 

(a) 


0.024377 

(b) 

Radial velocity 

7319 km s“^ 

(a) 


7308 km s"! 

(b) 

Distance 

104 Mpc 


Spatial scale 

504 pc arcsec”^ 


Diameter 

50 arcsec 25 kpc 

(b) 

R 25 

26.75T3.25 arcsec 

(c) 

Absolute rotation velocity 

lOOTll km s“^ 

(b) 

(at radii>22 kpc) 



Star formation rate 

2.19 M© yr-i 

(d) 


2.10T0.06 M© yr-i 

(b) 

Stellar mass 

1.86x10^0 M© 

(d) 

Hi mass 

2.3xlOio M© 

(e) 

(a)|Hickson et al.|(|1992|), (b) 

this work, (c)|de Vaucouleurs et al.| 

(|1991|), (d)|Bitsakis et al.|(|2014|), (e) |Borthakur et al. 

|(|2010|). 


Figure 1. Red band image of HCG 91a, c and d from DSS-2. The 
purple rectangles denote the footprint of our WiFeS observations 
of HCG 91c. The location of the sky field was chosen to avoid 
any of the faint and extended tidal structures stemming from 
HCG 91a, including the faint tidal arm extending (on-sky) from 
HCG 91a to HCG 91c. 


multi-phase, multi-scale view of HCG 91c to uuderstaud the 
preseuce of abuudauce aud velocity offsets iu some of the 
H H regious detected with WiFeS iuside this galaxy. 

The article is structured as follows. We first discuss our 
differeut datasets (iucludiug their acquisitiou, reductiou aud 
processiug) iu Sectiou We theu describe the group-wide 
structure aud kiuematic of the Hi gas iu HCG 91 iu Sec¬ 
tiou aud zoom-iu ou HCG 91c iu Sectiou]^ We discuss 
possible couuectious betweeu group-wide mechauisms aud 
specific star formiug regious iuside HCG 91c iu Sectiou 
aud summarise our couclusious iu Sectiou [H 

Throughout this paper aud uuless uoted other¬ 
wise, wheu we refer to specific emissiou Hues we 
meau [N ii] = [N ii] A6583, [S ii] = [S ii] A6717+A6731, 
[O II] = [O ii] A3727+A3729 aud [O iii] = [O iii] A5007. We 
adopt t he ma ximum likelihood cosmology from |Komatsu| 
Ho=70A km s”^ Mpc“\ Qa=0.73 aud 


et al. 


Qm=0.27, resultiug iu a distauce to HCG 91c of 104 Mpc. 


2 OBSERVATIONAL DATASETS 

2.1 WiFeS 

2.1.1 Observations 


We observed HCG 91c with WiFeS, the Wide-Field lutegral 
Field Spectrograph ( Dopita et al.||2007| 2010) mouuted ou 
the 2.3m telescope (Mathewsou et al. 2013) of the Australiau 
Natioual Uuiversity at Sidiug Spriug Observatory iu North- 
eru New South Wales, Australia. The footpriut of our ob- 
servatiou is showu iu Figure overlaid ou a red-baud image 
from the Second Digitized Sky Survey (DSS-2Q We com- 
biued two iudividual WiFeS poiutiugs to coustruct a 38x50 


arcsec^ mosaic iu order to better cover the spatial exteut of 
HCG 91c. 

The data was acquired over two uights ou 2012 August 
15 aud 16, with a seeiug of 1.2-1.5 arcsec. Each iudividual 
scieuce field of 25x38 arcsec^ was observed 4x 1400s, result¬ 
iug iu 5600s ou-source per field. A blauk patch of sky to the 
South-West of HCG Ola was observed for a total iutegra- 
tiou time of 700s betweeu every pair of cousecutive scieuce 
frames. 

WiFeS is a dual-beam iutegral field spectrograph with 
iudepeudeut red aud blue chauuels. We used the B3000 grat- 
iug for the blue arm aud the R7000 gratiug for the red arm, 
iu coujuuctiou with the RT560 dichroic. This setup allows 
us to have a complete waveleugth coverage from 3800A to 
7000A, as well as a high spectral resolutiou observatiou of 
the Ha emissiou hue to perform a detailed kiuematic aualy- 
sis. It is importaut to uote here that with this iustrumeutal 
settiug, the fiual spectra associated with each spectral pixel 
(a.k.a spaxel) iu the mosaic has a dual spectral resolutiou 
of R=3000 from 3800A to 5560A aud R=7000 from 5560A 
to 7000 A, correspoudiug to a full width at half maximum 
(FWHM) for emissiou hues of ~100 km s“^ aud ~43 km 
s“^, respectively. 


2.1.2 Data reduction 

Each scieuce exposure was iudividuahy bias subtracted, 
flat fielded, sky subtracted, waveleugth aud flux calibrated, 
as well as atmospheric refractiou corrected with pywifes 
(vO.5.6), the uew official python data reductiou pipeliue 
for WiFeS observatious (Childress et al. |2014a|b| | . This uew 
pipeliue eutirely replaces the uow obsolete IRAI0 pipeliue 


^ The DSS-2 data was obtained from the European Southern 
Observatory Online Digitized Sky Survey Server. 


^ IRAF is distributed by the National Optical Astronomy Obser¬ 
vatories, which are operated by the Association of Universities for 
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( Dopita et al.|20T0 | by providing several key improvements: 
scriptability, reproducibility, advanced wavelength calibra¬ 
tion technique using an accurate model of the instrument’s 
optics, multicore processing and compatibility with the first 
& second generations of CCD detectors in the WiFeS cam¬ 
eras. The pipeline is described in detail in [Childress et al.| 
(2014b), to which we refer the reader for more details about 


further details. We use the Gonzalez Delgado et al. (2005) 


the different functions. 

The reduced eight red science frames and eight blue 
science frames are median combined in two mosaics (red 
and blue) using a custom made PYTHON script. Because the 
wavelength sampling was chosen to be the same during the 
data reduction process, mosaicking only requires a shift of 
the data cubes in the spatial x- and y- directions. Given the 
seeing conditions and the WiFeS spaxel size of 1x1 arcsec^, 
we restrict ourselves to integer spatial shifts. The final mo¬ 
saic contains 38x50=1900 spaxels. Once the red and blue 
mosaics are constructed, we correct them for Galactic ex¬ 
tinction using the Schlafly & Finkbeiner ( |2011| ) recalibra¬ 
tion of the Schlegel et al.j ( |1998| extinction map based on 
dust emission measured by COBE/DIRBE and IRAS/ISSA. 
We assume for HCG 91c a Galactic extinction Aa/=0.052, 
following the NASA/IPAG Extragalactic Database (NED). 
The recalibration assumes a Eitzpatrick (1999) reddening 


law with Rv^S.l and a different source spectrum than that 
used by [Schlegel et al. ( 1998[ ). 

We note that the far blue end of a WiEeS spectrum be¬ 
low 4000A (with the B3000 grating) is a known problematic 
region. The core reason, aside of a drop of the overall in¬ 
strument transmission, is that the available flat-held lamps 
(at the epoch of our observations) have almost no flux at 
these wavelengths. It is therefore virtually impossible to ac¬ 
curately flat-held the data below 4000A with the B3000 grat¬ 
ing. PYWIFES combines lamp flat-held images with twilight 
sky flat-held images which mitigate this issue partially, but 
not entirely. In the present case, this flat-fielding issue di¬ 
rectly affects the [O II] flux, to which we And a correction 
factor of 1.5 ought to be applied. This correction factor is 
only an average for the entire held of view, so that we de¬ 
cided not to rely on the [O II] line flux in our spaxel-based 
analysis of HGG 91c. 


2.1.3 Spectral fitting 

We fit the strong emission lines for all 1900 spaxels in our 
mosaic after fitting and removing the underlying stellar con¬ 
tinuum in the spectra. This task is performed automatically 
using the custom build LZIFU vO.3.1 IDL routine, written by 

Ho 


I-Ting Ho at the University of Hawaii (Ho et al. 2014 


et ah, in prep.). LZIFU is a modified version of the UHSPEC- 
FIT IDL routine to which D.S.N. Rupke, J.A. Rich and H.J. 
Zahid all contributed. The code is a wrapper around PPXF 
( Gappellari Emsellem[[2004 ), used to perform the contin¬ 
uum fitting, and mpfit (Markwardt 2009| , used to fit the 
emission lines with different gaussian components. 

The LZIFU routine first performs (via ppxf) the contin¬ 
uum fitting by looking at spectral regions free of emission 
lines. We refer the reader to the PPXF documentation for 


Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. See http://iraf.noao.edu 


set of stellar templates with the Geneva isochrones to con¬ 
struct the best-matched stellar continuum for each spaxel. 
We note that we tested using the stellar templates based 
on the Padova isochrones instead, but found no significant 
difference in the fit results. 

LZIFU offers the possibility to fit one or multiple gaus¬ 
sian components to multiple emission lines simultaneously. 
In the present case, a visual inspection of the data revealed 
a very narrow structure of the different emission lines with 
no evidence of multiple peaks in any spaxel. For consistency, 
we performed two distinct fits of the emission lines - with 
one, and subsequently with two gaussian components. We 
detected the presence of a broadened base around the Ha 
emission line in the ~100 brightest spaxels, and confirmed 
the detection with a statistical f-test (e.g. Westmoquette[ 


et al. 2007). We performed this test blindly for the [Nil] 


and Ha emission lines only for all spaxels, and with a null- 
hypothesis rejection probability of 0.01. However, because 
a) the broadened base is visually only marginally detected 
around the Ha line in ~ 10 out of 1081 spaxels with signal- 
to-noise (S/N) larger than 5, b) it is not detected around 
any other emission line, and c) beam-smearing effects may 
be responsible for some of the broad component, we choose 
to enforce a 1-gaussian component fit for every spaxel in our 
mosaic of HGG 91c. 

We note here that nowhere in our observations do we de¬ 
tect clear multiple peaks in the emission line profile of Ha as 


reported by Amram et al. (2003) in their Fabry-Perot obser¬ 
vations of HGG 91c. Although the spectral sampling (0.35A) 
and spectral resolution (R=9375 at Ha) of their data are 
slightly better than ours (0.44A and R=7000), WiFeS (of 
which the instrumental design minimizes scattered light and 
other instrumental artefacts) should have clearly detected 
the distinct velocity components offset by ~100 km s“^ ob¬ 
served by Amram et al.[ ( [2003| . A detailed comparison of 
the WiFeS line profile with the velocity structures of [Amram| 
et al. (2003) is presented in Appendix [A| With no evidence 


for multiple components and/or significant asymmetries m 
the line profiles in our WiFeS observations, we are led to 
conclude that the gas in HGG 91c is associated to a single 
kinematic structure at all locations. 

The Ha flux map of HGG 91c observed by WiFeS is 
shown in Figure in units of 10~^^ erg s“^ cm“^. The final 
emission line widths quoted in this article are corrected for 
the instrumental resolution. 

Throughout this article, we analyze our WiFeS dataset 
on a per-spaxel basis to keep a high spatial resolution even 
at large distances from the galaxy center. To ensure trust¬ 
worthy results, we adopt varying S/N cuts throughout our 
analysis. For example, we require S/N(Ha;H/3)^5 to per¬ 
form the extragalactic reddening correction. The specific 
S/N cuts will be discussed individually for each case in the 
next Sections. For clarity, we also include the associated S/N 
cut and number of spaxel in all relevant diagrams through¬ 
out this article. 

A series of spaxels with S/N(Ha)^5 and S/N(H/3)<5 
are detected to the North of HGG 91c at the mean position 
[4;25]. This is the most distant star forming region (from the 
galaxy center) detected with WiFeS, and as such it holds 
unique clues regarding the state of the ionized gas at large 
radii. To ensure a reliable extragalactic correction, we have 
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Figure 2. Ha flux map of HCG 91c, before correcting the extra- 
galactic reddening. The core of the galaxy and the spiral structure 
is traced by enhanced Ha emission. Three star forming regions 
of interest “Cl”, “C2” and “C3” are labelled accordingly. The 
region “Rl” is also indicated, along with the individual spaxels it 
comprises. The x and y axis are in arcsec (= spaxels), centered on 
the brightest spaxel in Ha. At the distance of HCG 91c, 1 arcsec 
503 pc. 


summed all of the contiguous spaxels with S/N(Hq;)^ 5 into 
a combined spectrum, which we refer to as the “Rl” region. 
The individual Rl spaxels are marked with empty white 
squares in Figure We also indicate the location of three 
star forming regions with anomalous abundance and kine¬ 
matics offsets, the “Cl”, “C2” and “C3” regions. 

Extragalactic reddening 

We correct our WiFeS observations for extragalactic redden¬ 
ing using the line ratio. HCG 91c does not contain 

any active galactic nuclei (AGN) which could increase the 
intrinsic line ratio value Ra/ 3 , and we therefore assume that 
R(xp — 2.86 a s in a case B reco mbination (rie — 100 cm“^; 


T=10^ K, see Osterbrock|l989 ). We follow the methodology 

described in details in the A 
and use the extinction curve 
with Ry=4.5. The choice of 

ppendix A of Vogt et al. (2013), 

^ from Fischera & Dopita 

(2005) 

this extinction curve, very sim- 

ilar to that of Calzetti et al. 

(2000), is motivated by t] 

he fact 
in the 

sistent 

that Wijesinghe et al. (2011 

) showed that for galaxies 

GAMA survey ([Driver et al. 

2009), it provides self-con 


estimates of the star formation rates computed from the UV, 
Ha and [O III]. The resulting V-band extinction Ay in mag¬ 
nitudes is shown in Figure The extragalactic reddening 
(and correction) is only calculated (and applied) for the 530 
spaxels with S/N(Ha;H/3)^5, corresponding to a line inten¬ 
sity detection threshold of ~2xl0“^^ erg s~^ cm“^ A“^. 

The largest extinction is present towards the core of 
the galaxy with Ay > 2. Further out, the extinction drops 
to Ay?^l, with localized increases coherent with the spiral 


Figure 3. Ay extinction, calculated from the Ha/H^ emission 
line flux ratio. The center of the galaxy is the most affected by ex¬ 
tragalactic reddening, with Ay > 2. Enhanced extinction regions 
are also associated with the spiral structure detected in Figure 
The uncertainty in the Ay value is of the order of 0.6-0.8 for all 
spaxels. 


arms visible in Figure Unlike most star forming regions in 
the spiral arms, we detect no enhancement of the extinction 
associated with the Cl, C2 and C3 regions at the 1-sigma 
level (given an uncertainty of 0.6-0.8 in the Ay value). 


2.2 VLA 


The Hi distribution in HCG 91 was observed with the Very 
Large Array {VLA) under the program AV0285 (P.I.: L. 
Verdes-Montenegro) on 2005 October 5. The observations 
were carried-out in the VLA D-configuration with 3.25hr of 
on-source integration and dual polarization using two inter¬ 
mediate frequency (IF). The resulting total bandwidth of the 
datacube is 6.2 MHz (~1200 km s“^) with a resolution of 
97 kHz. The correlator was setup to cover the entire HI spec¬ 
tral linewidth with a velocity resolution of ~21.6 km s“^. 
The data was calibrated following the standard VLA cali¬ 
bration and imaging procedures using NRAO’s Astronom¬ 
ical Image Processing Software (AIPS). The final absolute 
uncertainty in the flux density scaling is ~10 per cent. The 
synthesized beam of the reduced datacube is 74” x47”, with 


a position angle of 4.83° West-of-North (Borthakur et al 


2010 ). 


2.3 Pan-STARRS 


HCG 91c was observed at multiple epochs with the Pan- 
STARR^l (PSl) telescope (Kaiser et al. 2002 Kaiser 2004 


Kaiser et al.|2010| [Hodapp et al.|200'4 ) on Haleakala (Maui) 


^ Panoramic Survey Telescope And Rapid Response System 
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as part of the PSl Stt survey (Chambers et ah, in prepa¬ 
ration). This survey mapped (several times) the entire sky 
visible from Hawaii ((5>-30°) in five broad filters (gpi:~400- 
550 nm; rpi:~550-700 nm; ipi: ^690-820 nm; zpi:~820-920 
nm; ypi :~920-1000 nm, see Tonry et al.|2012 ) between 2010 
and 2014. The Pan-STARRS images presented in Section [4T] 
were obtained from the online Postage Stamp server. These 
images were automatically processed by the PSl Image pro¬ 
cessing Pipeline (IPP, see e.g. Magnier||2007 Magnier et al. 


2013 ). Specifically, the Pan-STARRS images presented in 
this article (see Figure]^ have been obtained from the Stv 
survey processing version 2 (PV2) release. They have a pixel 
size of 0.25 arcsec, and are the result of the combination 
of 20 single epoch exposures obtained under natural seeing 
conditions (~1 arcsec). 


THE GROUP-WIDE HI DISTRIBUTION IN 
HCG 91 


Hi is an excellent tracer of past and ongoing gravitational 
interactions. For compact groups, HI can also be used to es¬ 


timate the evolutionary stage of the group overall (Verdes- 
Montenegro et al.|2001 ). The Hi map for HCG 91 obtained 
with the VLA is shown in Figure The left and right pan¬ 
els show Position-Position and Position-Velocity projections 
of the data cube, respectively. Iso-intensity surfaces are fit¬ 
ted to the distribution of HI emission in the 3-dimensional 


datacube using the MAYAVI module in python (Ramachan- 


dran Sv Varoquaux 2011), before being projected onto the 


respective 2-dimensional planes. Each iso-intensity level is 
shown with an 80 per cent transparency, with exception to 
the inner-most one with 0 per cent transparency. The loca¬ 
tions of the HCG Ola, b, c and d galaxies are marked with 
black crosses. The location of the Cl, C2 and C3 star form¬ 
ing regions in HCG 91c is indicated with small golden cubes 
and black crosses. HI structures of interest (see below) are 
marked with red-crosses for rapid identification. A large HI 
tail associated with HCG Ola is traced with a red rod. Fig- 
ure [ 4 ] is also inte ractive, following a concept described by 
Barnes & Fluke (2008). By clicking on it in Adobe Acro¬ 


bat Reader v9.0 or above, it is possible to load an interac¬ 
tive view of the (X;Y;v) data cube. The interactive model 
is also accessible via a supplementary HTML file compati¬ 
ble with most mainstream web browsers Q When inspecting 
this interactive 3-D map, one should remember that it is not 
“fully” spatial, but rather an X-Y-v volume. Presenting the 
VLA 3-D data cube as 3-D model is similar to the example 


provided by Kent (2013), although our respective 3-D model 
creation methods are different (i.e. we do not use blender). 
The total Hi mass associated with HCG 91 is 2.3xlO^°Mo 
( Borthakur et al.|2010 ). 

For ease of visualisation, we show in Figurethe differ¬ 
ent symbols marking the position of the galaxies, HI clumps 
of interest and tidal arm introduced in Figure]^ on top of 
a DSS-2 red band image, along with our WiFeS observation 
fields. The full extent of the different HI structures is traced 


At the date of publication, the interactive HTMLdocument is 
compatible with firefox, chrome, safari and internet ex¬ 
plorer. We refer the reader to the X3DOM documentation for 
an up-to-date compatibility list: http://www.x3dom.org/ 
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Figure 5. Red band image of HCG 91a,b,c and d from DSS-2. 
The image spans the same area as the Position-Position velocity 
diagram of Figure(left), of which the same markers are shown 
for comparison. The purple rectangles denote the footprint of our 
WiFeS observations of HCG 91c. 
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Figure 6. Red band image of HCG 91a,b,c and d from DSS-2. 
The image spans the same area as shown in Figure ^ (left) and 
Iso-intensity contours trace the Hi content of the group, and 
are color-coded as a function of their radial velocity, ranging from 
6822.0 km s“^ (dark blue) to 7405.5 km s“^ (dark red), in steps 
of '^21.6 km s“^. The contours are extracted from the different 
velocity slices of the VLA datacube and trace an intensity of 2.5 
mjy/beam. 


with iso-contours (at a level of 2.5 mJy/beam and colored 
as a function of the gas velocity) in Figure]^ 

Globally, the different sub-structures of the HI gas dis¬ 
tribution in the group can be associated with the different 
individual galaxies. Specifically: 

• HCG Ola: Spatially, this galaxy is coherent with the 
large HI structure to the South of the group. The large tidal 
feature detected in HI is also spatially coherent with a faint 
optical counterpart distinguishable in the red band DSS-2 
image of the area (see Figure]^. Kinematically however, we 
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Figure 4. Position-Position (left) and Position-Velocity (right) projection of the Hi gas content of HCG 91, as seen by the VLA. The 
gas density is traced by a series of iso-intensity semi-transparent surfaces (set at 1.3, 2.5, 3.5 and 6.0 mJy/beam, respectively). The 
location of the different group members is marked with black crosses, and labelled accordingly above either panels. Two Hi clumps 
associated with the main gas reserve around HCG 91c are marked with red crosses for unambiguous identification. The large tidal tail 
originating to the South-East of HCG 91a is traced by a red rod. The Cl, C2 and C3 star forming regions in HCG 91c are marked 
with golden cubes and black crosses (the golden cubes are 5 arcsec x 5 arcsec x 5 km s“^ in size). The corresponding VLA beam 
size (74”x47”, see [Borthakur et al.|[^10| > is shown to the bottom left of the left panel. An int e ractiv e v e rsion of this Figur e can b e 
acc e ss e d e ith e r with Adob e Acrobat R e ad e r v.9.0 or abov e by clicking on this Figur e , or via most mainstr e am w e b brows e rs by op e ning 

th e corr e sponding HTML fil e includ e d as suppl e m e ntary mat e rial to this articl e . Th e 3 - D cub e can b e fr ee ly rotat e d and/or zoom e d in 

and out. In th e int e ractiv e cub e , th e blu e axis is th e v e locity dim e nsion, th e r e d axis is th e R.A. dir e ction, and th e gr ee n axis th e D e c. 

dir e ction. A s e ri e s of pre-set vi e ws allow th e vi e w e r to e asily r e produc e proj e ctions of int e r e st or r e mov e th e low e st int e nsity iso - surfac e 

for a cl e ar e r vi e w of e ach galaxy. Th e gold e n cub e s marking th e location and spatial/kin e matic e xt e nt of th e Cl, C2 and C3 star forming 

r e gions ar e b e st s ee n in th e int e ractiv e count e rpart of this Figur e . NdlA: the interactive component of this Figure has been disabled in 
the arXiv version of the article (its compilation requires the mediaQ package, while arXiv currently only supports the obsolete movielh 
package). The interactive Figure will be accessible from the final article in MNRAS (both as an interactive PDF inside the article itself, 
and as a supplementary interactive HTML document). Until then, readers can access the interactive HTML version of this Figure at 
http:// WWW. mso. anu. edu. au/ ~ fvogt/website/misc. html . 


between the optical red- 


Hickson et al. 1992) and 


)pt 

T 


find a mismatch of ~220 km 
shift of the galaxy (6832 km s“^, se 
the mean redshift of the Hi structure. Although the optical 
redshift of HCG 91a associates the galaxy with a local Hi 
over-density in the Posit ion-Velocity diagram of Figure 
most of the Hi gas spatially-coherent with the galaxy (in¬ 
cluding the base of the tidal tail) is redshifted by up to ~400 
km s~^ with respect to the optical redshift. The Hi distri¬ 
bution observed with the VLA (both the spatial extent and 
the dynamic range of the different structures) is consistent 
with similar observations from Barnes Webster] ( |2001| us¬ 
ing the Australia Telescope Compact Array (ATCA), while 
the optical redshift measurement of HCG 91a from |Hickson| 
|et al. ( 1992| (z=0.022789) is consistent with the redshift 
measurement from the 6dF Galaxy Survey data release 3 
(z=0.022739, see [Jones et al.|20'Q4 2009) and with our own 
observations of the galaxy with WiFeS. Altogether, these 
observations confirm that the observed optical-radio redshift 
mismatch of HCG 91a is real. Here, we merely mention the 
existence of this redshift mismatch, and defer any further 


discussion to a later paper dedicated to our WiFeS observa¬ 
tions of HCG 91a. 

• HCG 91b: A spatially compact and kinematically ex¬ 
tended (~490 km s“^) Hi structure is associated with this 
galaxy. The large range of the Hi kinematics is consistent 
with the fact that HCG 91b is seen almost edge-on on. 

• HCG 91c: A rotating Hi disk with a velocity range 
of ~200 km s“^ (^7200-7400 km s“^) is associated both 
spatially and kinematically with the optical counterpart of 
HCG 91c. We find a small kinematic offset of ~10 km s~^ be¬ 
tween the optical redshift of HCG 91c measured by j Hickson] 


et al. (1992) and the mean velocity of the HI structure. The 


largely undisturbed morphology of the Hi distribution is 
suggesting the presence of only minimal tidal effects for this 
galaxy. To the North-West, two fainter Hi sub-structures 
(marked with red crosses in Figure]^ appear connected to 
the main gas reservoir of HCG 91c. They are also connected 
(less strongly) at the 1.3 mJy/beam level to the Hi struc¬ 
ture associated with HCG 91b. These two Hi clumps are 
located 116 arcsec ~ 58.4 kpc and 150 arcsec ~ 75.5 kpc 
from the center of HCG 91c. They have no visible optical 
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8 F.P.A. Vogt et al. 


counterpart in the DSS-2 red band image in Figure Spa- 
tially, HCG 91c is located ~15 kpc to the North-East of the 
large tidal arm originating from HCG 91a. The HI gas in the 
tidal arm is blueshifted by 150-250 km s“^ from the mean 
velocity of HCG 91c. 

• HCG 91d: This galaxy is not associated with any Hi 
structure kinematically. 


The large tidal feature originating from HCG 91a makes 
the HCG 91 compact group a Phase 2 group in the classifi¬ 
cation of Verdes-Montenegro et al. (2001), although some 
of the Hi gas is still clearly associated with the galax¬ 
ies HCG 91b and HCG 91c. The Hi reservoir associated 
with HCG 91c appears largely undisturbed from a kinematic 
point of view. The two HI gas clumps located to the North- 
West of HCG 91c may have resulted from tidal stripping, 
suggesting that HCG 91c is experiencing the first stages of 
tidal disruption via gravitational interaction. The HI bridge 
at the 1.3 mJy/beam level connecting the gas reservoir of 
HCG 91b and HCG 91c could be seen as evidence for an 
ongoing interaction between the HI envelopes of these two 
galaxies, although the exact bridge structure would require 
a higher spatial sampling to be clearly established. 


4 THE GALAXY HCG 91C 

Here, we describe the different characteristics of HCG 91c 
as seen by WiFeS and Pan-STARRS. We focus our analysis 
on the strong emission lines and the associated underlying 
physical characteristics of the ionized gas. We restrict our¬ 
selves to a description of the system, and postpone a global 
discussion of the different measurements until Section |5] 


4.1 Spatial distribution of the stellar light 


A color-composite image of HCG 91c, constructed from the 
Pan-STARRS ipi, rpi and gpi band images, is shown in 
Fig-IZl The same image is shown twice, with the intensity 
stretch adjusted to reveal the outer regions (left panel) and 
the inner spiral structure (right panel) of HCG 91c. The 
WiFeS mosaic footprint is traced by a dashed rectangle in 
the left panel, and the location of R 25 is traced by an el¬ 
lipse in the right panel. Specifically, R25=26.75±3.25 arcsec 
= 13.5±0.2 kpc with an ellipticity of 0.91, following NED 
and the Third Referenee Catalogue of Bright Galaxies (de 
Vaucouleurs et ar]|1991 ). We note that this value is consis¬ 
tent with the estimate of R 25 =28.45 arcsec and an ellipticity 
of 0.87 from the Surfaee Photometry Catalogue of the ESO- 
Uppsala Galaxies ( Lauberts V Valent ijn||l 989 ). 

HCG 91c harbors a regular, tightly wrapped spiral pat¬ 
tern with a bright nucleus. The outer limit of the stellar 
disk, traced by white iso-contours in the right panel of Fig.[^ 
(extracted from the gpi image) is mostly regular, with the 
brightness dropping the most sharply towards the SE. The 
spiral arms are extending (at least) 20 arcsec 10 kpc from 
the galaxy center, but their brightness decreases rapidly be¬ 
yond ~5 kpc. 

The Cl, C2 and C3 regions defined in Figure [^appear as 
three distinct compact features (each indicated by a white 
arrow in the left panel of Figure for clarity). The Pan- 
STARRS broad-band images primarily trace the stellar con¬ 
tent of HCG 91c, and one must use caution when comparing 


this view to that of the ionized gas provided by WiFeS. Nev¬ 
ertheless, with a diameter of ~ 6 Pan-STARRS pixels pe 1.5 
arcsec 750 pc across, it is very likely that the optical emis¬ 
sion associated with the Cl, C2 and C3 star forming regions 
is unresolved both in our WiFeS observations and in the 
Pan-STARRS image of HCG 91c. 


4.2 Emission line maps and line ratio diagrams 


We present the flux line maps for the principle optical emis¬ 
sion lines ([Oh], H/3, [OHI], Ho, [Nil], [Sll]) in the spec¬ 
trum of HCG 91c (before applying the extragalactic red¬ 
dening correction) in Figure]^ For comparison purposes, all 
the maps are shown with the same color stretch and cuts. In 
each panel, the white contours at 2.2xl0~^^ erg s~^ cm~^ 
trace approximatively the regions with S/N^5. 

The structure of the emission from the [N II] and [S II] 
lines is similar to Ho, although with an overall weaker in¬ 
tensity. There is a distinct lack of [O II] emission from the 
core of the galaxy, a region subject to a larger reddening (see 
Figure]^. The noise level is much higher in the [Oil] map, 
because the lines are located at the very blue end of the 
WiFeS spectra (in a region subject to flat-fielding problems, 
see Section 2.1.2). The spatial distribution of the [Olll] line 
is distinctly different from that of Ho. The lack of emission 
at the core of the galaxy is reminiscent of [O II], and is con¬ 
sistent both with a larger extinction and with a higher gas 
metallicity overall in the core of HCG 91c (see Section 4.3). 
Further out, the [O III] emission is the strongest in the Cl 
and C3 regions. By comparison, the Ho emission associated 
with the Cl and C3 regions is at a level comparable to other 
star forming regions within the spiral arms. 

The standard line ratio diagrams log [Nil]/Ho vs 
log [0HI]/H/3 and log [SH]/Ho vs log [OlU]/ll/3 for 


HCG 91c are shown in Figure Each 353 spaxels for which 
S/N(Ho;H/3)^5 and S/N([0 Hl];[N H];[S H])^2 are shown as 
individual squares, color-coded as a function of S/N([0 III]). 
The first S/N selection condition results from the extragalac¬ 
tic reddening correction described previously. A lower S/N 
cut for [O III], [N II] and [S II] is then acceptable given that 
both Ho and H^d are clearly detected and strongly constrain 
the velocity and velocity dispersion of all the different emis¬ 
sion lines. The color scheme choice reflects the fact that 
S/N([0IH])^ 2 is in all cases the limiting eondition. The 
grey error bars shown in each diagram indicate the 1-sigma 
error associated with the line ratio measurements, calcu¬ 
lated from the errors measured by LZIFU based on the vari¬ 
ance measurements propagated through the pywifes data 
reduction pipeline. 

Within the errors, the different spaxels form a well de¬ 
fined sequence in both line ratio diagrams, extending 1.5dex 
along the log [Olll]/ll/3 direction. The spaxel density along 
the sequence is continuous, except for a marked decrease 
at log [O Hl]/H/3?^0. This gap, best seen in the density his¬ 
togram in the right-hand side of the upper panel of Fig¬ 
ure is most certainly not a consequence of low S/N mea¬ 
surements. Indeed, spaxels on either side and across this 
gap in the distribution have strong oxygen emission with 
S/N([0 Hl]>6, and are therefore clearly detected. 

In Figure we show the log [O Hl]/H/3 line ratio map 
for all 353 spaxels visible in either panel of Figure Along 
with the R1 region, the spaxels associated with the Cl, C2 
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Figure 7. Color combination of Pan-STARRS ipi, rpi and gpi images (shown as the red, green and blue channels respectively) of 
HCG 91c. The color scheme is stretched to reveal the galaxy outskirts (left) or the inner region and spiral structure (right). The WiFeS 
observations footprint is marked with a dashed rectangle in the left panel, along with three arrows marking the location of the Cl, C 2 
and C3 regions. R 25 is traced by the dashed ellipse in the right panel. Also in the right panel, low level iso-contours constructed from a 
gaussian smoothed gpi band image highlight the full optical extent of HCG 91c on the sky at a distance of ~R 25 . 
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Figure 8. Maps for the different emission lines considered in this article, shown with a uniform color scheme and color stretch for 
comparison purposes. For each line, the white contours at 2.2x10“^^ erg s“^ cm“^ approximatively trace the regions with S/N^5. 


and C3 regions have log [O Hl]/H/3^0. They are forming 
the upper end of the star forming sequence visible in the 
line ratio diagrams. Altogether, Figures and indicate 
an abrupt change in the physical conditions of the ionized 
gas between the Cl, C2 and C3 regions and their immediate 
surroundings. 

4.3 Oxygen abundance and ionization parameters 

We now delve deeper into the physics of emission line ra¬ 
tios, and compute oxygen-abundances and ionization pa¬ 
rameters for each 353 spaxels with S/N(Ho;H/3)^5 and 


S/N([0 Hl];[N Il];[S H])^2. The concept and the method em¬ 
ployed here follow [Popita et al.| ( |2013a| . We use the pyqz 
vO.6.0 PYTHON module to compute, for each spaxel, the value 
of the oxygen abundance 12+log(0/H) and the ionization 
parameter log(q) (with q in cm s“^). The code relies on 
a series of grids of photoionization models computed with 
MAPPINGS IV - the most recent embodiment of the MAPPINGS 


code (Dopita et al.||1982| Binette et al.| 

11982 1985 ISuther- 

land & Dopita|1993| [Groves et al.|2004| 

1 [Allen et al.|2008 ) - 


which for a distinct set of line ratio diagrams allow to disen¬ 
tangle the values of 12+log(0/H) and log(q) over a large pa¬ 
rameter space (5 — 0.05 Zq and 6.5 ^ log(q) ^ 8.5). Among 
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Figure 9. Standard line ratio diagnostic diagrams 
log [Nii]/Hq; vs. log [Oiii]/H/5 (top) and log [SiiJ/Ha vs. 
log [Oiii]/H/3 (bottom). Each spaxel with S/N(H(a;H/3)^5 and 
S/N([0 lll];[Nll];[S ll])^2 are marked with squares colored as 
a function of S/N([Olll]). The position of the R1 integrated 
spectra is marked with a red circle and labelled accordingly. The 
1-sigma errors are shown in grey for each measurement. The 
full black lines are the maximum starburst diagnostic lines from 
[Kewley et al.| (|2001|>, the da shed line is the starburst diagnostic 
line from [Kauffmann et al.| (^03|>, while the Seyferts-LINERs 
separation from [Kewley et al. ( |2Q06| > is marked as a dash-dotted 
line in the bottom panel. The histogram in the upper panel shows 
the spaxel distribution as a function of log [Olll]/H/3, while the 
smooth distribution (traced by the black curve) was derived by 
performing a Kernel Density Estimation on the distribution of 
log [Oiii]/H/3 values. 


other updates, MAPPINGS IV allows for a non-Maxwellian 
distribution for the electron energies in the form of a Hi- 
distribution ( Nicholls et al. 12012 2013 Dopita et al.|2013b ). 
Here, we adopt At=20, but we note that this choice does not 
significantly affect our analysis, given the limited influence 
of a distribution on the intensity of strong emission lines. 

We note that for this work we have upgraded PYQZ from 
vO.4.0 publicly released by Dopita et al. (2013b) to vO.6.0. In 
this latest version, PYQZ can fully propagate errors in the line 
flux measurements through to the estimation of the abun¬ 
dance and ionization parameters (but still relies on the same 
MAPPINGS IV models of Dopita et al. [20131^ . The new PYQZ 
vO.6.0 will be publicly released in the near future along with 
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Figure 10. Map of log [Oiii]/H/3 for spaxels with S/N(Ho;H/3)^5 
and S/N([0 iii];[N ii];[S ii])^2. Most spaxels with a positive ratio 
are located in the Cl, C2 and C3 regions. The black line traces 
the iso-contour corresponding to log [O lll]/H/5=0. 


updated mappings iv models (Sutherland et ah, in prep.), 
and is available on-demand until then. For completeness, we 
describe in Appendix [B] how observational errors are propa¬ 
gated in this new version of PYQZ and how the final 1-sigma 
uncertainty level on the 12+log(0/H) and log(q) values is 
computed, namely via the propagation of the full probability 
density function associated with each emission line measure¬ 
ment. 

In total, eight diagnostic grids are available in PYQZ 
to compute the abundance and ionization parameter of a 
given spectrum. In principle, each of these grids can provide 
an estimate of 12+log(0/H) and log(q), which can then be 
combined to provide a global estimate. In practice, and for 
the present case, we only used the diagnostics not involving 
[O II] because these emission lines are subject to flat-fielding 
issues (see Section [2.1.2| . In addition, an inherent twist in 
the model grid over the region of interest for HCG 91c ren¬ 
ders the diagram log [N II] /Ha vs log [O IH]/Hy5 unusable in 
the present case. Ultimately, we are left with just two suit¬ 
able diagnostic grids: log [NH]/[SH] vs log [0HI]/H/3 and 
log [NH]/[SH] vs log [0HI]/[SH]. 

The resulting 12+log(0/H) and log(q) maps for 
HCG 91c are shown in Figure m As mentioned previ¬ 
ously, the overall decrease of the [O Hi] flux in the inner 
regions of HCG 91c is associated with increasing oxygen 
abundances. The peak oxygen abundance at the galaxy cen¬ 
ter is 12+log(0/H)?^9.25. By comparison, the oxygen abun¬ 
dance of the R1 region is found to be 12+log(O/H)?^8.80. 
The oxygen abundance of the Cl, C2 and C3 regions is of the 
order of 8.7-8.9. The ionization parameter map is uniform 
throughout the entire galaxy (7.10 < log(q) < 7.35). The Cl 
region is a notable and the only significant departure from 
the mean with log(q)~7.65. 

In Figure we construct the abundance gradient of 
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Figure 11. Gas phase oxygen abundance (top) and ionization 
parameter (bottom) measured with the PYQZ python module for 
all spaxels with S/N(Hq;;H/ 3)^5 and S/N([0 iii];[Nii];[S ii])^2. In 
the top panel, the dashed line traces the major axis of HCG 91c 
with a P.A. of 40 degrees West-of-North (see Section [4^. 


HCG 91c as a function of the deprojected radius. Every 
353 spaxels for which we derived an oxygen abundance are 
shown individually. The vertical error bar associated to each 
measurement corresponds to the error computed by PYQZ 
vO.6.0. To deproject the position of each spaxel in the disk 
of HCG 91c, we assume an ellipticity 5/a=0.8±0.08 (mea¬ 
sured manually from the rpi image of HCG 91c from Pan- 
STARRS), and a position angle P.A.=40° West-of-North 
(measured from the rotation map of HCG 91c, see Sec¬ 
tion This P.A. is in disagreement with the P.A. avail¬ 


able in NED based on near-IR images from 2MASS (Skrut- 
skie et aL||2006 ) of 5° West-of-North, but is in good agree¬ 


ment with the P.A. derived from optical images from the 


Surface Photometry Catalogue of the ESO-Uppsala Galax¬ 
ies ( Lauberts fc Valentijn||1989 ) of 46° West-of-North. The 
horizontal error bar associated with each measurement cor¬ 
respond to an estimated 10 per cent measurement error on 
the ellipticity. We assumed the center of the galaxy to be 
located at the position [0,0], coincident with the peak Ho 
emission, which given the seeing conditions is also consistent 
with the kinematic center of the galaxy, as described in the 
next Section. 

As illustrated in the top panel of Eigure |12| a linear 
fit to the full set of data points fails to properly match 
the outer regions of the disk. It is clear from Eigure m 
that our sampling of the outer regions of the disk (where 
12+log(0/H)^8.92) is poor, and largely composed of the 
Cl, C2 and C3 regions. Hence, the global gradient fit is 
largely influenced by the spaxels in the inner regions of 
HCG 91c. In the middle panel of Eigure we use a lin¬ 
ear broken fit to better reproduce the trend in the oxy¬ 
gen abundance throughout HCG91c. The fit is broken at 
4.5 kpc from the galaxy center - the approximate radius at 
which the oxygen abundance gradient appears to steepen. 
An inner gradient is derived from spaxels in the Zone 1 (Zl; 
radius<4.5kpc, 12+log(0/H)>8.92), while an outer gradi¬ 
ent is derived from spaxels in the Zone 2 (Z2; radius^4.5kpc, 
12+log(0/H)>8.92), and forced to match the inner gradient 
at 4.5 kpc. The respective slopes of the different gradients 
are compiled in Table All oxygen abundance gradients 
were derived using the Orthogonal Distance Regression rou¬ 
tines inside the SCIPY module in PYTHOr^ 

A series of spaxels beyond 7 kpc from the galaxy center 
remain below the best fit outer gradient by ~0.15 dex. In the 
bottom panel of Eigure we show in black the spaxels lo¬ 
cated within 3 arcsec from the galaxy’s major axis visible in 
Eigure These spaxels trace an abrupt drop in the oxygen 
abundance of ~0.2 dex at 7 kpc from the galaxy center over 
a distance of 1 kpc = 2 WiEeS spaxels. This sharp oxygen 
abundance decrease is essentially unresolved (spatially) in 
our data set, so that its true nature (physical discontinuity 
in the oxygen abundance or sharp but continuous decrease) 
remains uncertain. It is however clear that the intensity of 
the oxygen abundance drop at 7 kpc is inconsistent with 
any linear gradient extrapolated from the inner regions of 
HCG 91c. 


We deliberately use only spaxels in the zone Z2 to de¬ 
rive the best-fit linear outer gradient, as we detect spaxels 
with 12+log(0/H)>8.92 for (almost) all azimuths beyond 
a radius of 4.5 kpc. Our sampling of less enriched gas is 
much less uniform, so that including all the spaxels with 
12+log(0/H)<8.92 in the outer gradient fit would wrongly 
bias it towards the Cl, C2 and C3 regions of which the 
abundances (as illustrated in the bottom panel of Eigure [T^ 
are clearly inconsistent with a linear gradient. Although the 
outer gradient is derived solely from the Z2 spaxels, we note 
that its slope nevertheless matches some of the outer regions 
of HCG 91c, including the R1 region. 

In the top panel of Eigure we color-code the dif¬ 
ferent spaxels as a function of their associated S/N([0HI]. 
Most spaxels lying in the Zone 4 (Z4; 12+log(0/H)<8.92, 


^ http://docs.scipy.org/doc/scipy/reference/odr.html, accessed 
on 2015 February 12 
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Figure 12. Oxygen abundance gradient in HCG 91c. Each 353 
spaxels for which 12+log(0/H) could be measured (see Figure 
are shown individually with their associated 1-sigma errors. The 
galaxy center, taken to be the peak of the Ho; emission line flux 
map, is at the position [0,0] (see Figure]^. A position angle P.A. 
= 40° and an ellipticity b/a=0.8 are assumed to deproject the 
radius of each spaxel within the disk of HCG 91c. The horizontal 
error bars correspond to a measurement error of 10 per cent of 
the ellipticity. The vertical error bars are the error computed by 
PYQZ vO.6.0. The linear best fit (and associated 1-sigma error) 
to all the spaxels is shown in the top panel. A broken linear fit, 
better matching the outer regions of the galaxy, is shown in the 
middle panel. In the bottom panel, the spaxels located within 
3 arcsec from the galaxy’s major axis (shown in Figure [TTJ are 
drawn in black with their associated errors in red. They reveal 
the abrupt drop of ~0.2 dex in oxygen abundance at ~7.5 kpc to 
the North-West of the galaxy center: the location of the Cl, C2 
and C3 star forming regions. 


Table 2. Characteristics of the different oxygen abundance gra¬ 
dients presented in Figure [T^ 


Gradient 

Fitted region(s) 

Slope 

[dex kpc“^] 

Zero point 

global fit 

all 

-0.029T0.001 

9.242T0.005 

inner fit 

Z1 

-0.020T0.003 

9.219T0.007 

outer fit 

Z2 

-0.038T0.002 

9.300T0.011 


12+log(0/H) more than 0.075 dex below the best fit outer 
gradient based on the Z2 spaxels) are clearly detected, and 
their errors (shown in dark red for clarity) place them 1-2 
sigma below our best-fit outer gradient. Their offset in oxy¬ 
gen abundance is best seen in the Figure’s bottom panel, 
showing the distribution of oxygen abundances for all spax¬ 
els in the zones Z2 and Z3+Z4. While some spaxels below 
12+log(0/H)=8.92 are consistent with the linear gradient 
estimate based on the Z2 spaxels out to 13 kpc (i.e. spaxels 
in the zone Z3), most lie 0.15 - 0.3 dex lower (in the zone 
Z4). 


Most spaxels in the zone Z4 belong to the C1,C2 and 
C3 regions, as illustrated in Figure It can be noted that 
as the Cl, C2 and C3 regions are located along or near the 
P.A. of HCG 91c, their deprojected radii are not subject to 
large uncertainties. The Z3 spaxels are on the other hand 
all located in the outer regions of the spiral arm extending 
Northward from the left-hand side of HCG 91c. Especially, 
the lack of S/N in the [-10;15] area (see Figure 11) appears 
consistent with the lack of detection of Z3 spaxels at a radius 
of 10-11 kpc. 


4.4 Kinematics of the ionized gas 


In Figure 15 we show the ionized gas velocity map (cen¬ 
ter) and velocity dispersion (right) of HCG 91c. In the left 
panel, we show the Ho intensity map overlaid with the iso¬ 
velocity contours. We use the freely available PYTHON rou¬ 
tine FIT_KINEMATIC_PA from M. CappellarJ^ to measure the 
P.A. of HCG 91c from its gas kinematics: P.A. = 40 zb 4 
degrees West-of-North. This routine (originally written in 
idl) was used extensively by Cappellari et al.| (2007) and 
Krajnovic et al. (2011), and is described in Appendix C 


of Krajnovic et al. (2006). The P.A.=40° (black & white 
dashed line) direction and the galaxy center (white dot) used 
to construct the abundance gradient of HCG 91c shown in 
Figure are shown in the left and middle panel of Figure p^ 
for completeness. 

The kinematic signature of the gas in HCG 91c is con¬ 
sistent with regular rotation. The quoted velocitie s assume 
a rest frame velocity of 7319 km s“^ (z=0.024414; Hickson 
et ah! 1992). We find an overall asymmetry between the red- 


shifted and blueshifted sides that could be reconciled if the 
rest-frame velocity was reduced by 11 km s“^ to 7308 km 


® http://www-astro.physics.ox.ac.uk/~mxc/software/, accessed 
on 2015 February 12 
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Figure 15. Velocity (center) and velocity dispersion (right) maps of HCG 91c, for all spaxels with S/N(Hq;)^5. Iso- velocity contours 
are shown in the middle panel, and are also overlaid over the Ha flux map of HCG 91c is the left panel. In the left and middle panel, 
the center of the galaxy (traced by the peak Ha line flux) is marked with a white circle, and the black & white dashed line traces the 
principal axis of HCG 91c (P.A.=40° West-of-North). The dark-red & white dashed line is rotated by 12° counter-clockwise from the 
galaxy’s major axis, and is designed to pass through the C1-C2-C3 complex of star forming regions. The Position-Velocity diagrams in 
Figure [T g] are extracted along these two directions. For the iso-velocity contours, the thick-full lines are spaced by 30 km s“^, the thin-full 
lines are spaced by 15 km s“^ and the dashed lines are spaced by 5 km s“^. 


s“ (as measured by fit_kinematic_pa). This offset would 
also help recoucile the optical redshift of HCG 91c with the 
HI kiuematic of the galaxy measured by the VLA (see Sec- 
tiou|^. 

The velocity dispersiou of the iouized gas is through¬ 
out the eutire disc of HCG 91c raugiug from ~20 km s“^ 
to 40 km s“^. Accouutiug for the thermal broadeuiug of the 
Hues (ctth ~ 13(T/10^)°‘^ km s“^ for the Hydrogeu hues, 
see Osterbrock||1989] ), this correspouds to au iutriusic ve¬ 
locity dispersiou rauge of ~ 15-38 km s“^. The largest ve¬ 
locity dispersious are fouud towards the galaxy ceuter, aud 
are most certaiuly iuflueuced by beam smeariug, giveu the 
seeiug couditious duriug our observatious (1.2-1.5 arcsec). 
No siguificaut iucrease of the velocity dispersiou is detected 
towards the Cl, C2 aud C3 clumps. 

We detect some deviatious from a regular rotatiou sig- 
uature at aud arouud the positious [5,15] to [12,12] - the 
locatiou of the Cl, C2 aud C3 regious. Specifically, the Cl, 
C2 aud C3 regious are redshifted from (i.e. are laggiug be- 
hiud) the regular rotatiou of the galaxy by 5-10 km s“^. We 
recall that the error associated with the spectral calibratiou 
of t he red Wi FeS datacube is of the o rder of O.OSA (see Sec- 
Childress et al.||2Q14b ) or 2.3 km at Ha, 


tiou 


2.1.2 


aud 


so that a 5-10 km s“^ offset correspouds to a 2-3 sigma de- 
tectiou. The fact that the kiuematic distortiou is exteuded 
aud tracks the bright spaxels distributiou of the Cl, C2 aud 
C3 regious suggest that the distortiou is real, although a 
greater spectral resolutiou aud fiuer spatial sampliug would 
be beueficial to coufirm its existeuce. 


To further quautify the distortious of the velocity map 
of HCG 91c associated with the Cl, C2 aud C3 star formiug 
regious, we coustruct two Positiou-Velocity (PV) diagrams. 
We choose the project iou axis to be orieuted a) aloug the 
galaxy’s P.A. aud b) aloug au axis rotated by 12° couuter- 
clockwise from the galaxy’s major axis. lu both cases, we 
set the zero-poiut at the galaxy ceuter aud add 11 km s“^ 
to the measured velocities to correct the global kiuematic 
asymmetry meutioued previously. The value of 12° bears uo 
special siguificauce other thau eusuriug that the secoud axis 
passes through the C1-C2-C3 complex of star formiug re¬ 
gious. The resultiug PV diagrams are showu iu Figure 
Every spaxel withiu 3 arcsec aud 1 arcsec (respectively) from 
the two projectiou axis aud with S/N(Ha)^5 for which the 
Hue velocity cau be measured accurately is showu as au iu- 
dividual red square (redshifted side of the galaxy) or a blue 
diamoud (blueshifted side of the galaxy), as a fuuctiou of its 
(projected) distauce aloug the PV axis. 

Aloug the major axis of the galaxy, after a liuear 
iucrease out to ~2 arcsec, the slope of rotatiou veloc¬ 
ity decreases aud remaius cousisteut betweeu the red- aud 
blueshifted sides out to ~12 arcsec. While the redshifted gas 
velocity flatteus out at ~75 km s“^ beyoud 12 arcsec, the 
blueshifted gas iucreases up to 80zb3 km s“^ beyoud 22 arc¬ 
sec. Assumiug au ellipticity for HCG 91c of 0.8zb0.08, we fiud 
the absolute rotatiou velocity of the gas at a radii greater 
thau 11 kpc to be lOOzbll km s“^. 

A uoticeable feature of the top diagram is the kiuematic 
behaviour of the blueshifted gas from 12 to 20 arcsec. As the 
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Figure 13. Top: same as the middle panel of Figure |12| but 
with the individual spaxels colored as a function of S/N([Olll]). 
The diagram is divided in four gradient zones (labelled Z1 to 
Z4) based on the overall trend morphology. For clarity, errors 
for spaxels with S/N([0 lll])^9 are shown with thicker dark-red 
lines. Bottom: histogram and associated smooth distribution of 
the oxygen abundance offset from the best-fit outer-gradient, for 
spaxels in the Z2 zone (black) and in the Z3+Z4 zones (red). 


magnified inset diagram in Figure [^illustrates, two velocity 
“branches” separated by 5-10 km s“^ at the 1-2 sigma level 
exist in this range. Especially, the “bottom” branch corre¬ 
sponds to spaxels in the C3 region, and is the signature of 
the lag of the C3 star forming region mentioned previously. 
The velocity lag of the Cl, C2 and C3 complex of star form¬ 
ing regions is best seen in the bottom panel of Figure [Tb] in 
which the blue velocity curve decrease by ~5-10 km s^ at 
16 arcsec from the galaxy center (i.e. at the location of the 
C2 star forming region). 

Spiral galaxies in compact groups have been observed 
to host a wide range of rotation signatures with some¬ 


times large asymmetries and/or perturbations (see e.g. Ru¬ 


bin et al. 1991). As compact groups favor strong gravita¬ 


tional interactions between galaxies, the existence of small 
perturbations in the velocity field of HCG 91c (located in¬ 
side a compact group) is in itself not surprising. Of interest 
however is the fact that the star forming regions with local¬ 
ized lower oxygen abundances are associated with localized 
kinematic anomalies. This abundanceF4kinematic connec- 


S/N(Hcx:H13)>5 

S/N([Olll];[Nll]:[Sll])>2 
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Figure 14. Spatial distribution of the spaxels in the different 
gradient zones Z1 to Z4 defined in Figure [T^ 


tion indicates that both the gas composition and motion in 
the Cl, C2 and C3 star forming regions are inconsistent with 
their immediate surroundings. 


4.5 Star formation rate 


We can convert the Ho line flux of each spaxel in our WiFeS 
mosaic (see Figure]^ into a star formation rate (SFRna) 
following the r ecipe of [Murphy et al. (2011) derived using 
STARBURST99 ( Lcitherer et al.|1999| ) and assuming a Kroupa 
(2001) initial mass function (IMF): 


= 5.37 X 10-42 ^ 

ergs“^ 


( 1 ) 


where Lua is the total Ho luminosity, derived from our mea¬ 
sured de-reddened fluxes per spaxel and given the assumed 
distance to HCG 91c of 104 Mpc. This recipe results in a 
SFR 68 per cent of what would be derived using the lKennu] 


cutt 


_ ^(1998) formula, largely because of differences in the as¬ 
sumed IMF characteristics ( Calzetti et al.|2007||Kennicutt 
Evans|20T2 ). The resulting SFR map of HCG 91c is shown in 
Figure 17 where only the 530 spaxels with S/N(Ho; 


and the R1 region that were corrected for extragalactic red¬ 
dening are shown. We note that our spatial resolution of 
0.5 kpc arcsec”^ is similar to the scale of the star forming 


complexes studied by Murphy et al. (2011), so that issues 


associated with using a “mean” Ha-to-SFR conversion fac¬ 


tor on small scales does not apply in our case (see Murphy 
et al.||2011 Kennicutt fc Evans||2012 ). 


As expected from the Ha emission line map, star for¬ 
mation activity in HCG 91c is most intense in the central re¬ 
gion of the galaxy, with additional hot spots along the spiral 
arms. The total SFRna for HCG 91c (summed from all the 
spaxel with S/N(Ha; H/3)^5 and the R1 region) is 2.10±0.06 
Mq yr“^. This value is in excellent agreement with the es- 
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Figure 17. Star formation rate derived from the Ho: flux, for all 
spaxel corrected for extragalactic reddening. The galaxy center 
is the most active region, and the spiral structure is also visible. 
Localized enhancements of the star formation rate are detected 
in the Cl, C2 and C3 regions. 


Table 3. Main integrated characteristics of the Cl, C2 and C3 
star forming regions 


Figure 16. Top: Position-Velocity diagram of the velocity map 
shown in Figure extracted along the axis passing through the 
galaxy center and with P.A. = 40° West-of-North (shown with 
a black & white dashed line in Figure p^ . Each spaxel within 3 
arcsec of the axis and with S/N(Hq;)^ 5 is shown as a red square 
(redshifted side of HCG 91c) or a blue diamond (blueshifted side 
of HCG 91c). A correction of 11 km s“^ was added to the ve¬ 
locity of each spaxel to correct the global kinematic asymmetry 
detected in the velocity held. Bottom: idem, but for spaxels with 
S/N(Hq;)^ 5 located within 1 arcsec from an axis rotated by 12° 
counter-clockwise from the galaxy’s major axis (traced with a 
dark-red & white dashed line in Figure p^ . 


timated 2.19 M© yr“^ from Bitsakis et al. (2014) based on 
a full SEE[^ fitting of the integrated light of HCG 91c. 

Assuming a total stellar mass = 1.86x10^° M© for 
HCG 91c ( [Bitsakis et al. 2014), we obtain a specific star 


formation rate sSFR = 1.13zb0.05x 10 yr This is com¬ 


parable to the lowest sSFR measured by [Tzanavaris et al.] 
(2010) in 21 star forming HCG galaxies from their UV and 


IR colors, but well within the distribution of sSFR reported 
by Plauchu-Frayn et al. (2012) for ~50 late-type galaxies 
(SBc and later) in compact groups using full optical spec¬ 
trum htting with stellar population synthesis models. In 
fact, HCG 91c is extremely consistent with the SFR vs M* 


^ Spectral Energy Distribution 


Region 

Cl 

C2 

C3 

< AR.A. > 

[ arcsec ] 

3.2 

7.0 

10.5 

< ADec. > 

[ arcsec ] 

13.5 

13.5 

15.3 

V gas (line of sight) 

[km s-i ] 

7245A4 

7236A3 

7230A4 

Cgas 

[km s-i ] 

21±2 

20±2 

20±2 

F(Ho) 

[IQ-ie gj.g g-i cjn-2j 

64.0A8.9 

23.8A6.0 

49.6A9.0 

SFRna 

[10-2 Mq yr-1] 

4.4A0.6 

1.7A0.4 

3.4A0.6 

Lna 

[1039 gj.g g-lj 

2.1A0.3 

0.8A0.2 

1.6A0.3 


and sSFR vs M* relations constructed from the SDSS and 
GAMA datasets for z < 0.1 star forming galaxies by [Lara-] 
Lopez et al. (2013). For reference, we list in Tablej^the dif¬ 


ferent characteristics of the Cl, C2 and C3 regions, whilst 
their integrated spectra are presented in Figure p^ 
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Figure 18. Raw integrated spectra (black) and LZIFU fit (red) of t' 
present in each panel are labelled on top of each column for clarity. 



Cl, C2 and C3 star forming regions. The dominant emission lines 


5 DISCUSSION: THE PECULIAR STAR 
FORMING REGIONS OF HCG 91C 


By and large, HCG 91c appears as a rather unremarkable 
star forming spiral galaxy. It hosts regular star formation 
rates throughout its disc comparable to other similar galax¬ 
ies in the field, and an overall regular rotation (as traced 
by the ionized gas). This galaxy also hosts a linear oxygen 
abundance gradient out to ~ 4.5 kpc, which then breaks to 
a steeper linear gradient. The Hi gas distribution around 
the galaxy traced by the VLA is largely undisturbed. The 
cold dust contours measured by Herschel/SPIRE at 250 /xm 
display a high degree of azimuthal symmetry ( [Bitsakis et al. 


2014). Mendes de Oliveira et al.| ( |2003| report that HCG 91c 


is an outlier with respect to the B-band Tully-Fisher rela¬ 
tion, in that its absolute B-band magnitude Mb is 2 mag 
brighter than expected from its rotation velocity, and con¬ 
cluded (also based on the two distinct kinematic components 
reported by Amram et al.||2003 ) that HCG 91c may be the 
result of a merger. We disagree with this picture after find¬ 
ing no evidence for the existence of multiple emission line 
component in our WiFeS dataset. Furthermore, the clear 
and regular spiral structure detected with WiFeS and Pan- 
STARRS strongly argues against the merger scenario, and so 
is the regular photometric profile of the disk (already noted 
by Mendes de Oliveira et al.||2Q03 ). 


Altogether, these characteristics suggest that HCG 91c 
has not (yet) strongly interacted gravitationally with the 
other galaxy members of HCG 91.The presence of faint, ex¬ 


tended HI gas to the North-West of the main HI reservoir 
(see Figure]^ suggests that HCG 91c may be caught in the 
very early stage of its interaction with the group. We note 
that the iso-intensity contours tracing the optical extent of 
the stellar disk of HCG 91c (see Figure reveal a sharp 
intensity drop to the South-East, but a smoother and more 
irregular boundary to the West & North-West, in the direc¬ 
tion of the kinematic asymmetries and HI clumps. 

Beside these large scale features, we found abundance 
and kinematic anomalies (at the 1-2 sigma level) at the 
location of the Cl, C2 and C3 star forming regions. Our 
emission line analysis reveals an abrupt change in the phys¬ 
ical condition of the ionized gas between these star forming 
regions and their immediate surroundings. They are com¬ 
paratively more metal poor by 0.15 dex, and the Cl re¬ 
gion hosts the largest value of the ionization parameter in 
the disk of HCG 91c with log(q)?^7.65. Such large values of 
log(q) can be explained by a very young age of the associated 
star forming region (<0.5 Myr), a l ow pressure environ ment 
(log(P//c)?^4-5 cm“^ K), or both (Dopita et al. 2006). We 
also find evidence that the gas in the Cl, C2 and C3 regions 

may be lagging behind the rotation of the disk by 5-10 km 
^-1 


The lack of S/N is our WiFeS observations beyond 7 
kpc from the galaxy center hinders us from assessing the de¬ 
tailed behaviour of the oxygen abundance gradient at these 
large radii. It is clear that the oxygen abundance of the 
Cl, C2 and C3 regions are clearly inconsistent with a linear 
gradient extrapolated from the inner regions of the galaxy 
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(see Figure 12). On the other hand, the oxygen abundances 
along the spiral arm extending Northward from the top-left 
of HCG 91c are found to be consistent with the extrapo¬ 
lated linear oxygen abundance gradient. This suggests that 
the oxygen abundance drop associated with the Cl, C2 and 
C3 star forming regions is localized to these locations, rather 
than indicative of a uniform drop in the abundance gradient 
for all azimuths. Determining whether the Cl, C2 and C3 
star forming regions are unique, or whether more star form¬ 
ing regions in the outer regions of HCG 91c display similar 
(localized) abrupt decrease in oxygen abundances will re¬ 
quire deeper follow-up observations (see Section 5.3). 


5.1 Origin of the fuel for star formation 

How can we explain these localized, anomalous Cl, C2 and 
C3 star forming regions ? Their comparatively lower oxygen 
abundance (i.e. the presence of less enriched gas) leads us 
to formulate three possible scenarios to explain their origin: 

(a) the accretion of a smaller, near-by satellite system, 

(b) gas inflow onto the disk of HCG 91c from the galaxy 
halo, and 

(c) the infall and collapse of pre-existing neutral gas clouds 
at the disk-halo interface of HCG 91c. 

For clarity, we discuss each of these scenarios separately. 


5.1.1 Satellite aeeretion 


The accretion of a near-by dwarf galaxy could possibly ex¬ 
plain the lower metallicity of the Cl, C2 and C3 regions. 
Lower metallicity dwarf galaxies have been detected around 
the Milky Way and M31, and the presence of large scale stel¬ 
lar streams implies the ongoing accretion and tidal stretch¬ 


ing of some of these systems (e.g. Belokurov et al. 2006 
Richardson et ar]|2011 ). The apparent on-sky alignment of 


the Cl, C2 and C3 regions (see Figure]^ is certainly sug¬ 
gestive of a physteal connection. Such a physical connection 
could then be explained by a tidally disrupted structure. 
However, the gas kinematic implies difficulties for this sce¬ 
nario. The observed lag of 5-10 km s“^ is small compared 
to typical infall velocities of streams in the halo of M31 and 
around the Milky Way. Although not impossible, this sce¬ 
nario would require a very specific set of circumstances be¬ 
tween the galaxy’s orientation, the accretion trajectory, and 
our point-of-view to result in the small velocity lag of the 
Cl, C2 and C3 star forming regions. 


5.1.2 Gas inflow from the galaxy halo 

The spatial and kinematic structure of the HI gas associated 
with HCG 91c certainly support the idea that this galaxy’s 
halo (still) contains a gravitationally-bound gas reserve. Un¬ 
der these circumstances, the lower oxygen abundance of the 
Cl, C2, and C3 regions could be explained if less enriched 
and unstructured material from the halo is in-falling onto 
the disk, and fuelling localized star formation activity. The 
small velocity lag and the low velocity dispersion of the gas 
in the Cl, C2 and C3 regions would suggest a slow infall 
velocity. Yet, the lack of any detectable enhancement of the 
velocity dispersion in the Cl, C2 and C3 regions compared 


to other star forming regions in the disk of HCG 91c remains 
puzzling. The exact triggering mechanism for the gas infall 
also remains an open question, but the presence of near¬ 
by galaxies and their gravitational fields appear as a likely 
source of perturbations. 


5.1.3 Infalling and eollapsing gas elouds at the disk-halo 
interfaee 

In the halo of the Milky Way, there exists a complex mix 
of gas in different phases, including molecular and neutral 
hydrogen ( Putman et al.|[2012 ). Part of this gas is located 
in a series of compact clouds with a wide range of kinemat¬ 
ics ( Saul et ^|2Q12 ). Some of these clouds have measured 
velocities largely inconsistent with galaxy rotation, and are 
usually refereed to as High Veloeity Clouds (HVCs, see e.g. 
Wakker van Woerden|[l997 ). By comparison. Intermedi¬ 
ate Veloeity Clouds (IVCs, |Wakker|2(30T ) also have velocities 
inconsistent with galactic rotation, but less so than HVCs. 
In fact, there exists a wide range of Hi cloud characteris¬ 
tics in the Milky Way’s halo, some of which are found to be 
co-rotating with the disk (see [Kalberla V Kerp 2009 and 
references therein). There also exists diffuse Hi gas in the 
halo of the Milky Way which displays a vertical lag in its 
rotation velocit y with respect to the disk o f the order of 15 
km s“^ kpc“^ ( Marasco V Fraternali||2011 ). 

Resolving the structure of the HI gas in other galaxies 
is observationally challenging. In recent years, the HALO- 
GAS survey (Heald et al. 2011) has revealed the presence 


of an extended Hi disk around edge-on, near-by galaxies 
( Zschaechner et al.||2011 [Gentile et al.||2013| |. These 

observations complement older detection of gaseous halos. 


for example around NGC 891 by Oosterloo et al. (2007) 


Most of these detections revealed a vertical lag of the Hi 
rotation speed above the galaxy disks of the same order of 


magnitude than in the Milky Way (but see also Kamphuis 
et al.|2013| |. 

Given the star forming, spiral nature of HCG 91c, it 
does not appear unreasonable to assume that it possess a 
multi-phase, complex gaseous halo similar to that of the 
Milky Way. Especially, although we do not know its detailed 
structure, the VLA observations clearly reveal a large rotat¬ 
ing Hi reservoir. Under these circumstances, could the Cl, 
C2 and C3 star forming regions have resulted from the in¬ 
fall and subsequent collapse of pre-existing gas clouds at 
the disk-halo interface of HCG 91c ? The compact nature of 
these star forming regions could be naturally explained by 
a “localised cloud” origin. Their velocity lag of 5-10 km s“^ 
would suggest that these star forming regions are in fact not 
located within the disk of HCG 91c, but 300-700 pc above it, 
if one a ssumes a vertical rotational v elocity lag of 15 km s“^ 
kpc“^ ( Marasco V Fraternali|[20lT ). Of course, the face-on 
nature of HCG 91c makes it impossible to directly measure 
any vertical offset for these star forming regions. However, 
we note that the associated reddening is lower in the Cl, 
C2 and C3 regions (Ay ^ 0.9, see Figure than in any 
other star forming regions in the spiral arms of HCG 91c 
(Av ^ 1.2) at the 1-sigma level. A lower reddening value 
may be the result of the location of the Cl, C2 and C3 
regions above the main disk (and dust) of HCG 91c. 
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5.2 Star formation triggering mechanism 


Clearly, we cannot firmly rule out any of the above scenarios 
invoked to explain the origin of the anomalous star forming 
regions Cl, C2, and C3 in HCG 91c. At this stage, we favor 
the idea that these star forming regions originated in the 
infall and subsequent collapse of pre-existing gas clouds at 
the disk-halo interface. Indeed, collapsing pre-existing gas 
clouds at the disk-halo interface could naturally explain the 
different properties of the anomalous star forming regions 
(velocity lag, lower oxygen abundance, compactness). 


The precise mechanism able to trigger the collapse of 
neutral gas clouds at the disk-halo interface of HCG 91c 
(and the subsequent rapid formation of molecular gas to fuel 
star formation) remains undefined. In the Milky Way, molec¬ 


ular hydrogen was detected in several IVCs (Richter et al. 
20Q3| Wakker|[2QQ6 ), and it has been proposed that 


tral hydrogen is compressed into molecular hydrogen as gas 
clouds fall onto the galaxy disk and get compressed via ram 
pressure stripping ( Weifi et al.|1999| Gillmon Shull| 2006 


Rohser et al.|2014 ). The Cl, C2 and C3 star forming regions 


may have resulted from a boosted version of this mechanism, 
where the initial infall of neutral gas clouds onto the galaxy 
disk is being triggered by large scale gravitational pertur¬ 
bations from near-by galaxies in an harassment-like process 
( Moore et al.||1998 ). The existence of tidal perturbations to 
the North-West of HCG 91c (supported by the existence of a 
possible HI tail and the comparatively complex edge of the 
stellar disk to the North-West of HGG 91c) may have given 


rise to local tidal shears or compressive tides (Renaud et al. 


2008 2009), triggering the collapse of the clouds. Theoreti¬ 


cally, Renaud et al. (2014) observed the effect of compressive 
turbulence and how it can lead to starburst events in in¬ 
teracting galaxies. This mechanism is already active during 
the early phase of galaxy interactions, and could therefore 
be active in HGG 91c. In the same simulations, large scale 
gas flows across a galaxy’s disk only occur at later stages 
of the interaction (i.e. at the second closest approach). We 
note that the initial star formation triggered by compressive 
turbulence does not appear to require large velocity disper¬ 
sions (cr <20 km s“^), which would be consistent with our 
WiFeS observations. 


Alternatively, HGG 91c is located ~30 arcsec = 15 kpc 
(on sky) to the North-East of the extended tidal tail of 
HGG 91a (see Figure]^. HGG 91c’s Hi envelope also sug¬ 
gest that the galaxy may have previously interacted with 
HGG 91b (as indicated by the presence of a possible Hi 
bridge between the two galaxies). Ram-pressure from either 
interaction may have shocked and compressed HGG 91c’s 
halo, leading to the collapse of the Gl, G2 and G3 star form¬ 
ing regions. Kinematically, HGG 91c and the tidal tail from 
HGG 91a are offset by ~250 km s“^, but given their spa¬ 
tial location, it is possible that HGG 91c is currently inter¬ 
acting/colliding with this extended tidal structure. Such a 
collision is impossible to firmly rule out without X-ray ob¬ 
servations of the group which may reveal the presence of hot, 
shocked gas resulting from the interaction. In any case, the 
relatively undisturbed HI envelope of HGG 91c and the con¬ 
tinuous structure of the tidal tail stemming from HGG 91a 
detected by the VLA would suggest that their interaction is 
at an early stage. 


5.3 Constraining the exact nature of the Cl, C2 
and C3 regions 

Our WiFeS observations, combined to VLA and Pan- 
STARRS datasets, have revealed the anomalous oxygen 
abundance of three compact star forming regions in the disk 
of HCG 91c, and showed that this galaxy is only just be¬ 
ginning its interaction with the compact group HCG 91. 
From the ionized gas physical characteristics and kinemat¬ 
ics, we found that infalling (and subsequently collapsing) gas 
clouds at the disk-halo interface could naturally explain the 
anomalous nature of the Cl, C2 and C3 regions. However, 
one should keep in mind that our different pieces of evidence 
for the anomalous nature of the Cl, C2 and C3 star forming 
regions (velocity and abundance offset, lower extinction) are 
all at the 1-2 sigma level. 

We also lack critical pieces of information: for example, 
a deeper understanding of the state of the gas (i.e. pressure, 
density, temperature) in these specific star forming regions, 
as well as a better characterization of the underlying stel¬ 
lar population. Additional insight on the stellar population 
associated with the Cl, C2 and C3 regions is especially crit¬ 
ical to test their possible origins. For example, collapsing 
halo gas clouds with masses ranging from 10^-10^ M© (see 
Putman et al. 2012 and references therein) do not offer a 
sustainable source of fuel and cannot host a significant pop¬ 
ulation of old stars, of which the detection would be a very 
strong argument against this scenario. Unfortunately, the 
S/N in the continuum of our observations does not allow us 
to extract meaningful information on the underlying stellar 
population in HCG 91c (see Figure p^. 

The inherent compact aspect of the Cl, C2 and C3 re¬ 
gions would also benefit from a finer spatial sampling to 
reveal their precise structural extent and possible physical 
connections. Finally, the disk of HCG 91c extends beyond 
the field-of-view of our WiFeS observations, and may con¬ 
tain additional anomalous star forming regions. To address 
these questions, we have been awarded 4.0 hr of Science 
Verification time on MUSE ( [Bacon et al.||2010 ), the new in¬ 
tegral field spectrograph on the Yepun telescope (unit 4 of 
the VLT) at Paranal in Chile (P.I.: F.P.A. Vogt, P.Id.: 60.A- 
9317[A]), which will be the subject of a separate article. 

We have shown that overall, HCG 91c is a largely undis¬ 
turbed (yet) star forming spiral galaxy. Hence, HCG 91c of¬ 
fers us a window on the early stage of galaxy interactions in 
a compact group, and possibly on the early phase of galaxy 
pre-processing in these environments. The presence of lo¬ 
calized star forming regions with lower oxygen abundances 
indicate that lower-metallicity gas is being brought from the 
outer regions of the disk or the halo to the inner regions of 
the galaxy. Hence, this gas may eventually contribute to the 
flattening of the oxygen abundance gradient in the system. 
Strong interactions of galaxies in pairs have been observed to 
flatten the metallicity gradient in these system through large 
scale gas flows ( Kewley et al.|20To Rupke et al.|2010 ). Our 
WiFeS observations of HCG 91c suggest that prior to large 
scale gas flows induced by strong gravitational perturbation, 
some gas mixing can be occurring at the level of individual 
star forming regions as a result of galaxy harassment and 
longer range gravitational interactions. 

Evaluating the ubiquity (or not) of collapsing gas clouds 
at the disk-halo interface of galaxies as an early consequence 
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of galaxy harassment will require additional observations of 
a statistical sample of galaxies. Existing or upcoming IFU 
surveys such as SAMI or MaNGA may provide such a statis¬ 
tical sample probing a wide range of environment densities. 
However, with a respective spectral resolution of R=4500 
and R=2000 and a spatial sampling ^ 1 kpc, finding sub-kpc 
star forming regions with ~0.15 dex oxygen abundances off¬ 
sets will certainly prove challenging for these surveys. Rather 
than “direct” detections, the SAMI and MaNGA surveys 
may instead provide “tentative” detections of anomalous 
star forming regions in largely unperturbed star forming disk 
galaxies. These objects would form an ideal sample (span¬ 
ning a wide range of environments) for dedicated follow¬ 
up observations at higher spectral and spatial resolution. In 
fact, we expect that large IFS surveys will provide an es¬ 
sential database allowing the identification of largely undis¬ 
turbed galaxies (similar to HGG 91c) in the early stage of 
their interaction with a compact group or cluster. As yet 
mostly unaffected by their environments, these systems can 
shine a unique light on the early stages of complex gravita¬ 
tional interactions and the associated consequences. 


6 SUMMARY 

In this article, we presented the discovery of three compact 
star forming regions with oxygen abundance and kinematic 
anomalies (at the 1-2 sigma level) in the otherwise unre¬ 
markable star forming spiral galaxy HGG 91c. From the 
analysis of the different strong optical emission lines de¬ 
tected with WiFeS, we found these anomalous star forming 
regions a) to be comparatively more metal poor by 0.15 dex 
with respect to their surrounding, and as expected from the 
overall metallicity gradient present in the inner region of 
HGG 91c, b) to kinematically lag behind the disk rotation 
by 5-10 km s“^, and c) for one of them, to be associated with 
the highest value of the ionisation parameter in the entire 
galaxy (log(q)?^7.95). 

To understand the origin of these peculiar-star forming 
regions, we combined our WiFeS data set with broad-band 
images of HGG 91c from Pan-STARRS, and with VLA ob¬ 
servations of the group-wide HI distribution of HGG 91. 
These datasets reveal that HGG 91c is still largely undis¬ 
turbed, but most likely experiencing the very first stage of its 
gravitational interactions with other galaxies inside HGG 91 
(and possibly with the large tidal tail from HGG 91a). 

Under these circumstances, we discussed three possible 
scenarios to explain the origin of the anomalous star form¬ 
ing regions detected with WiFeS: accretion of a satellite, gas 
inflow from the halo, and collapsing pre-existing gas clouds 
at the disk-halo interface. We found that the latter scenario 
could naturally explain all of the observed characteristics of 
the anomalous star forming regions (lower metallicity, ve¬ 
locity lag, compactness, lower reddening). By comparison, 
the satellite accretion and gas inflows scenarios are harder 
to reconcile with the observed gas kinematics, but cannot 
be firmly ruled out at this stage. 

We discussed possible mechanisms able to trigger star 
formation in these originally stable, pre-existing gas clouds, 
in the form of tidal shears, long-range gravitational pertur¬ 
bations or harassment from the other group members. The¬ 
oretically, these scenarios are consistent with the recent sim¬ 


ulations of Renaud et al. (2014) suggesting that compressive 
turbulence is responsible for enhanced star formation activ¬ 
ity in the very early stages of galaxy interactions, while large 
scale gas flows within a galaxy disk feeding central starbursts 
occur later on at the second closest approach. 

HGG 91c may be offering a direct window on the early 
phase of galaxy interactions in Gompact Group environ¬ 
ments, and possibly on one of the early stage of galaxy 
evolution. The existence of localized regions of lower metal¬ 
licity gas suggest that gas mixing (leading to a flattening 
of the overall abundance gradient) can be occurring on the 
scales of individual star forming regions before the onset of 
large inflows following stronger gravitational effects. In the 
era of large scale IFS surveys such as GALIFA, SAMI and 
MaNGA, HGG 91c act as a reminder that mechanisms asso¬ 
ciated with galaxy evolution may first be impactful on sub- 
kpc scale and display a discreet kinematic signature (Av?^10 
km s“^ in the present case). 

Dedicated follow-up MUSE observations will provide 
us with a sharper view of HGG 91c and its peculiar star 
forming regions. This dataset ought to let us better char¬ 
acterise the physical conditions of the ionized gas through¬ 
out HGG 91c, refine our detections of lower abundances, 
kinematic offsets and reduced extragalactic reddening asso¬ 
ciated with the Gl, G2 and G3 star forming regions, as well 
as provide important information regarding their underly¬ 
ing stellar population. The MUSE observations will be key 
to test our current favoured scenario for explaining the lower 
oxygen abundances of the Gl, G2 and G3 regions, and let 
us identify whether additional star forming regions in the 
outer regions of HGG 91c (undetected with WiFeS) display 
similar gaseous abundance anomalies. 
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APPENDIX A: COMPARISON OF THE GAS 
KINEMATICS WITH AMRAM ET AL. (2003) 


Amram et al. (2003) presented Fabry-Perot observations of 


the Ha emission line in the galaxy HGG 91c, and reported 
the existence of large asymmetries and multiple distinct 
components in the line profile (see their Figure 20). Their 
data was acquired with the GIGALE Fabry-Perot mounted 
on the ESO 3.6m telescope at La Silla on 1995 August 21-24 
with a seeing of ~1 arcsec. The spectral resolution of their 
dataset was R=9375 at Ha, with 24 scanning steps and a 
sampling step of 0.35A or ~16 km s“^, and a spatial pixel 
size of 0.91x0.91 arcsec^. 

In Figure M we reproduce the Figure 20 of |Amram| 


et al. (2003) showing the Ha line profile for the central 
15 xl5 arcsec^ regions of HGG 91c, but for our WiFeS data. 
Modulo minor differences (i.e. each panel corresponds to 1 x 1 
arcsec^ instead of 0.91x0.91 arcsec^), both Figures are di¬ 
rectly comparable. Especially, each panel in our Figure m 
covers the same spectral range as Figure 20 of|Amram et al. 
( [ 200 ^ . 

These authors reports secondary velocity components 
redshifted by ~100 km s“^ from the main line peak, and 
with dispersion extending over the entire spectral range 
shown in the different panels. It is clear from Figure [AT] that 
despite the slightly lower spectral resolution (R=7000) and 
spectral sampling (0.44A or ~19.5 km s“^) of our dataset, 
WiFeS would have detected complex line profiles such as 
those reported by [Amram et al. ( 2003| . Instead, the Ha 
line profiles observe d by WiFeS are narrow (with a 20-40 
km s“^, see Section 4.4) and single-peaked. We detect some 
small asymmetries in the inner most spaxels consistent with 
beam smearing. We note that these asymmetries are mostly 
blueshifted with respect to the main line peak, while the 
asymmetries reported by Amram et al. (2003) are largely 
redshifted. 

As we find no evidence for multiple kinematics compo¬ 
nents in our WiFeS observations of HGG 91c - and especially 
find no evidence for the complex line profiles reported by 


Amram et al. (2003) - we are lead to conclude that the gas 


in this galaxy can be well described by a single kinematic 
structure at all locations. 


APPENDIX B: ERROR PROPAGATION IN 
PYQZ VO.6.0 


In the original version published by Dopita et al. (2013b), 
the PYQZ vO.4.0 module was not set up to account for errors 
in the flux measurements of emission lines. Instead, errors on 
the estimation of log(q) and 12+log(0/H) were estimated 
from the standard deviation between the different estimates 
stemming from the different line ratio diagnostic grids cho¬ 
sen by the user. Indeed, for flux measurements at the ~5% 
level, the mismatch between the estimates from different di¬ 
agnostics typically dominates the uncertainty of the joint 
final estimate. For larger observational errors on the emis¬ 
sion line fluxes, PYQZ vO.4.0 therefore under-estimates the 
true error associated with the final values of 12+log(0/H) 
and log(q). Here, we have upgraded PYQZ to address this 
shortcoming. The updated version of PYQZ will be publicly 
released and described in details in the near future along 
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Figure Al. Reproduction of Figure 20 of Amram et al. (2003), but for our WiFeS observations: the Ho: line profile for all the spaxels 
within a box of 15x15 arcsec^ centered on the galaxy’s core (i.e. spaxel [0;0] in our reference frame). Each panel contains 19 wavelengths 
bins ranging from 6719.56A to 6727.dSA. Each panel is normalized to the emission line peak. 


with new MAPPINGS IV models. Here, we briefly describe 
for completeness how errors on the final 12+log(0/H) and 
log(q) estimates in our analysis of HCG 91c are computed. 


The underlying idea is to propagate the full probability 
density function associated with each line flux measurement 
to the 12+log(0/H) vs log(q) plane. To that end, for each 
set of line fluxes (and errors) provided by the user (i.e. for 
each individual spaxel, in the case of our analysis), PYQZ 
vO.6.0 generates n random set of line fluxes from the errors 
provided by the user, assuming a normal and uncorrelated 
error distribution for each line. PYQZ subsequently derives an 
estimate of 12+log(0/H) and log(q) for all n set of random 
line fluxes, and for every line ratio diagnostic grid selected 
by the user, resulting in n x m estimates of 12+log(0/H) 
and log(q) (where m is the number of line ratio grids se¬ 
lected by the user). Effectively, these n x m estimates rep¬ 
resent the discretized two-dimensional joint probability dis¬ 
tribution function of the oxygen abundance and ionization 
parameter. The full joint probability density function can 
the be reconstructed by performing a two-dimensional ker¬ 
nel density estimation. This frequentist approach to error es¬ 
timation differs from the Bayes formalism adopted by |Blan^ 


et al. (2014) in the izi code written in IDL with a similar 


purpose to PYQZ. 


It is often practical to derive from the joint probability 
density function a best-estimate and associated uncertainty. 
To derive the best estimate for 12+log(0/H) and log(q), 
PYQZ normalizes the reconstructed joint probability func¬ 
tion to its peak level, and computes the 61%-level contour 
(corresponding the 1-sigma level for a normal distribution). 
The mean location of this contour is then chosen as the best 
estimate of the 12+log(0/H) and log(q), while the full ex¬ 


tent of the contour is indicative of the uncertainty associated 
with these estimates. 

The process is illustrated for two representative high 
and low S/N spaxels in our WiFeS observations of HCG 91c 
in Figures B1 and |B2| respectively. For reasons detailed 
in Section 4.3 these examples rely on only two line ra¬ 
tio diagnostic grids: log [NH]/[SH] vs log [OlH]/H/3 and 
log [NH]/[SH] vs log [O HI]/[S II]. The location of the in¬ 
dividual estimates stemming from either diagnostic grids in 
the 12+log(0/H) vs log(q) are marked with white squares. 
The mean of these positions, corresponding to the PYQZ 
vO.4.0 estimate of 12+log(0/H) and log(q), is marked with 
a white star. For both spaxels, 1000 random set of line 
fluxes were generated by PYQZ, resulting in 2000 estimates 
of 12+log(0/H) and log(q) (1000 for either line ratio grid), 
shown as individual red dots. 

The 61%-of-peak-level contour of the underlying recon¬ 
structed joint probability density function is shown with 
an orange line. With 1000 random sets of line fluxes, we 
find that the location of the peak of the joint probability 
density function (marked with an orange diamond symbol) 
varies within the 1-sigma contours when PYQZ is run multi¬ 
ple times. However, the location of the 1-sigma contour at 
61% of the peak level remains consistent when PYQZ is run 
multiple times. This fact motivates our choice to base the 
final estimates of 12+log(0/H) and log(q) on the mean lo¬ 
cation of the 1-sigma contour level, rather than the exact 
peak location. 

Different PYTHON tools exists to perform a two- 
dimensional kernel density estimation. PYQZ can be set by 
the user to use two of these - either the GAUSSIAN_KDE rou- 
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tine in the SCIPY. STATS package or the kdemultivariate 
routine in the STATSMODEL package The former routine 
is 10-100 times faster but is not well suited for bi- or multi¬ 
modal distributions. In the case of PYQZ, the GAUSSIAN_KDE 
routine tends to over-smooth the joint probability function 
when two diagnostics are not directly consistent with one an¬ 
other (for example in the case of a flux measurement error). 
The option to use the GAUSSIAN_KDE routine is included in 
PYQZ vO.6.0 as a rapid alternative to test the code and ob¬ 
tain preliminary results. In practice, it is preferable to use 
the KDEMULTIVARIATE routine instead which provides more 
accurate results, as the bandwidth of the gaussian kernel 
can be set individually for the 12+log(0/H) and log(q) di¬ 
rections. All estimates of oxygen abundances and ionization 
parameters presented in this article were derived using the 
KDEMULTIVARIATE method. 


® http://docs.scipy.org/doc/scipy-0.14.0/reference/generated/ 
scipy.stats.gaussian_kde.html, accessed on 2014 December 13. 

® http://statsmodels.sourceforge.net/devel/generated/ 
statsmodels.nonparametric.kernel_density.KDEMultivariate.html, 
accessed on 2014 December 13 
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Figure Bl. Full (left) and magnified (right) oxygen abundance and ionization parameter plane accessible to PYQZ vO.6.0. The data 
corresponds to the spaxel [0;0] in our observations of HCG 91c. The white squares mark the two estimates associated with the diagnostic 
grids log [Nll]/[Sll] vs log [Olll]/H^ and log [Nll]/[Sll] vs log [O lll]/[S ll]. The red dots mark the location of 2x1000 randomly- 
reconstructed line fluxes from the original line fluxes and associated errors. The reconstructed 2-D joint probability density distribution 
(with a fixed resolution of 0.01 in 12+log(0/H) and 0.01 in log(q)) is also shown, with the orange contour tracing the 1-sigma level from 
the peak. The peak location is marked with an orange diamond, while the mean location of the 1-sigma contour and its extent in marked 
with the orange circle and associated error bars. The dashed box marks the full extent of the random set of individual estimates (the 
red dots). 



Figure B2. Same as Figure [ bT] but for the spaxel [9;14]. The lower S/N in the data results 
estimates. 


in a larger spread of the random set of 
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